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Abstract 

We present an algorithm for multi-hop routing and scheduhng of requests in wireless networks in 
the SINR model. The goal of our algorithm is to maximize the throughput or maximize the minimum 
ratio between the flow and the demand. 

Our algorithm partitions the links into buckets. Every bucket consists of a set of links that have 
nearly equivalent reception powers. We denote the number of nonempty buckets by a. Our algorithm 
obtains an approximation ratio of 0{a ■ log n), where n denotes the number of nodes. For the case 
of linear powers cr = 1, hence the approximation ratio of the algorithm is 0(log n). This is the first 
practical approximation algorithm for linear powers with an approximation ratio that depends only 
on n (and not on the max-to-min distance ratio). 

If the transmission power of each hnk is part of the input (and arbitrary), then a — 0(log F + 
log A), where F denotes the ratio of the max-to-min power, and A denotes the ratio of the max-to- 
min distance. Hence, the approximation ratio is 0(log n ■ (log F + log A)). 

Finally, we consider the case that the algorithm needs to assign powers to each link in a range 
[-Pmin , ^niax] ■ An cxtcnslon of the algorithm to this case achieves an approximation ratio of O [(log n+ 
loglogF)-(logF + logA)]. 

1 Introduction 

In this paper we deal with the problem of maximizing throughput in a wireless network. Throughput 
is a major performance criterion in many applications, including: file transfer and video streaming. It 
has been acknowledged that efficient utilization of network resources requires so called cross layered 
algorithms IILSS061 . This means that the algorithm deals with tasks that customarily belong to different 
layers of the network. These tasks include: routing, scheduhng, management of queues in the nodes, 
congestion control, and flow control. 

The problem we consider is formulated as follows. We are given a set V of n nodes in the plane. 
A link e is a pair (sg, re) of nodes with a power assignment P^. The node Se is the transmitter and 
the node is the receiver. In the SINR model, receives a signal from Sg with power Se = Pe/d^, 
where is the distance between Sg, and r^. and a is the path loss exponent. The network is given 
a set of requests {Ri}^^^. Each request is a 3-tuple Ri = (sj,tj,6j), where Sj € F is the source, 
ii £ V is the destination, and hi is the requested packet rate. The output is a multi-commodity flow 
/ = (/i, • • • , /fc) and an SiNR-schedule S = {Lt}JSQ that supports /. Each Lt is a subset of links that 
can transmit simultaneously (SiNR-feasible). The goal is to maximize the total flow |/| = Yli=i \fi\- 
We also consider a version that maximizes minj=i ^ \fi\/bi. Let A = djnux/dmm is the ratio between 
the maximum and minimum length of a link, and F = Pmax/^min the ratio between the maximum 
and minimum transmission power. For the case in which maxg^e' ;f7 = 0{1), the approximation ratio 
achieved by the algorithm is O(logn). For arbitrary powers and link lengths, the approximation ratio 
achieved by the algorithm is 0(log n • (log F + log A)). 
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Previous Work. Gupta and Kumar IGKOOH studied the capacity of wireless networks in the SINR- 
model and the graph model for random instances in a square. The SiNR-model for wireless networks was 
popularized in the algorithmic community by Moscibroda and Wattenhofer iMW061 . NP-Completeness 
for scheduling a set of links was proven by Goussevskaia IIGOW07 1. 

Algorithms for routing and scheduling in the SiNR-model can be categorized by four main criteria: 
maximum capacity with one round vs. scheduling, multi-hop vs. single-hop, throughput maximization 
vs. latency minimization, and the choice of transmitter powers. In the single-hop setting, routing is not 
an issue, and the focus is on scheduling. If the objective is latency minimization, then each request is 
treated as a task, and the goal is to minimize the makespan. 

The following problems are considered. (1) CAP-lSLOT: find a subset of maximum cardinality that 
is SINR-feasible. (2) LAT-lHOP: find a shortest SINR-schedule for a set of hnks. (3) L AT- PATHS: find a 
shortest SINR-schedule for a set of paths. (4) LAT-ROUTE: find a routing and a shortest SINR-schedule 
for a set of multi-hop requests. (5) THROUGHPUT- ROUTE: find a routing and maximum throughput 
SINR-schedule for a set of multi-hop requests. We briefly review some of the algorithmic results in this 
area published in the last three years. 

Chafekar et al. [C KM+O?! present an approximation algorithm for LAT-ROUTE. The approximation 



ratio is 0(log n • log A • log^ F). Fanghanel et al. MFKVIOI improved this result to 0(log A • log^ n). 
Goussevskaia et al. IIGWHW091 pointed out that log A can be Q{n), and presented the first approxima- 
tion algorithm whose approximation ratio is always nontrivial. In fact, the approximation ratio obtained 
by Goussevskaia et al. |,GWHW09J is O(logn) for the case lat-Ihop with uniform power transmis- 
sions. 

Halldorsson HHalOQI presented algorithms for LAT-lHOP with mean power assignments. He pre- 
sented an 0(logn log log A) -approximation and an O (log A) -online algorithm that uses mean power 
assignments with respect to OPT that can choose arbitrary power assignments (see also UTonlOII ). 

Halldorsson and Mitra BHMllall presented a constant approximation algorithm for CAP- 1 SLOT prob- 
lem with uniform, linear and mean power assignments. In addition, by using the mean power assignment, 
the algorithm obtains a 0(log n+log log A)-approximation with respect to arbitrary power assignments. 

Kesselheim and Vocking |KV10] give a distributed randomized algorithm for LAT-lHOP that ob- 
tains an 0(log^ n)-approximation using uniform and linear powers. Halldorson and Mitra BHMllbH 
improve the analysis to achieve an 0(log n)-approximation. 

Kesselheim [Kesll] presents approximation results in the SiNR-model: an 0(l)-approximation for 
CAP-lSLOT, an O (log n) -approximation for LAT-lHOP, an 0(log^ n) -approximation for LAT-PATHS 
and LAT-ROUTE. In liKesIll there is no limitation on power assignment imposed neither on the solution 
nor on the optimal solution. In practice, power assignments are limited, especially for mobile users with 
limited power supply. 

The most relevant work to our result is by Chafekar et al. lICKM+OSt who presented approximation 
algorithms for THROUGHPUT-ROUTE. They present the following results, an 0(log A)-approximation 
for uniform power assignment and linear power assignment, and an 0(log F • log A) for arbitrary power 
assignments. 

For linear powers. Wan et al. llWFJ+llI obtain a 0(log n) -approximation for THROUGHPUT-ROUTE. 
The algorithm is based on a reduction to the single-slot problem using the ellipsoid method. In liWanOQI . 
Wan writes that "this algorithm is of theoretical interest only, but practically quite infeasible." For the 
case that the algorithm assigns powers from a limited range. Wan et al. jWFJ"*"!!! achieve an 0(logn ■ 
log r)-approximation ratio. 



Our result. We present an algorithm for THROUGHPUT-ROUTE. Our algorithm partitions the links 
into buckets. Every bucket consists of a set of links that have nearly equivalent reception powers. We 
denote the number of nonempty buckets (also called the signal diversity of the links) by a. Our algorithm 
obtains an approximation ratio of 0(cr ■ log n), where n denotes the number of nodes. 



2 



For the case of linear power assignment the signal diversity is o" = 1, hence the approximation ratio 
of the algorithm is O(logn). This is the first practical approximation algorithm for linear powers that 
obtains an approximation ratio that depends only on n (and not on ratio of the max-to-min distance). 
This improves the 0(log A)-approximation of Chafekar et al. lICKM+OSt for linear power assignment. 
As pointed out in IIGWHW091 . log A can be Q{n). The linear power assignment model makes a lot 
of sense since it implies that, in absence of interferences, transmission powers are adjusted so that the 
reception powers are uniform. 

In the case of arbitrary given powers, the signal diversity is cr = 0(logr + log A). Hence, 
the approximation ratio is 0(logn • (logF + log A)). For arbitrary power assignments Chafekar et 
al. lICKM+OSi presented approximation algorithm that achieves approximation ratio of 0(log F • log A). 
In this case, the approximation ratio of our algorithm is not comparable with the algorithm presented by 
Chafekar et al. lICKM+OSt (i.e., in some cases it is smaller, in other cases it is larger). 

For the case of limited powers where the algorithm needs to assign powers between Pmin and Pmax> 
we give a 0[(log n + log log F) • (log F + log A)]-approximation algorithm. 

Our results apply both for maximizing the total throughput and for maximizing the minimum fraction 
of supplied demand. Other fairness criteria apply as well (see also liChaOQJ ). 



Techniques. Similarly to lICKM+OSt our algorithm is based on linear programming relaxation and 
greedy coloring. The linear programming relaxation determines the routing and the flow along each 
route. Greedy coloring induces a schedule in which, in every slot, every link is SiNR-feasible with 
respect to longer links in the same slot. 

We propose a new method of classifying the links. In lICKM+OSi IHal091 the links are classified by 
lengths and by transmitted powers. On the other hand, we classify the links by their received power. 

We present a new linear programming formulation for throughput maximization in the SiNR-model. 
This formulation uses novel symmetric interference constraints, for every link e, that bound the interfer- 
ence incurred by other links in the same bucket as well as the interference that e incurs to other links. 
We show that this formulation is a relaxation due to our link classification method. 

We then apply a greedy coloring procedure for rounding the LP solution. This method follows IIABL05[ 



CKM+08 , Wan09| and others (the greedy coloring is described in Section [63] ). 

The schedule induced by the greedy coloring is not SiNR-feasible. Hence, we propose a refinement 
technique that produces an SiNR-feasible schedule. We refine each color class using a bin packing 
procedure that is based on the symmetry of the interference coefficients in the LP. We believe this 
method is of independent interest since it mitigates the problem of bounding the interference created by 
shorter links. 



Organization. In Sec.|2]we present the definitions and notation. The throughput maximization prob- 
lem is defined in Sec.|3] In SeclH we present necessary conditions for SiNR-feasibility for links that are 
in the same bucket. The results in Sec. |4]are used for proving that the Unear programming formulation 
presented in Sec. |5]is indeed a relaxation of the throughput maximization problem. The algorithm for 
linear powers is presented in Sec. |6]and analyzed in Sec. |7] In Sec. [8] we extend the algorithm so that it 
handles arbitrary powers. In Sec. |9]we extend the algorithm so that it assigns hmited powers. 



2 Preliminaries 

We briefly review definitions used in the literature for algorithms in the SINR model (see l|HW09[ 
ICKM+081 ). 

We consider a wireless network that consists of a set V of n nodes in the plane. Each node is 
equipped with a transmitter and a receiver. We denote the distance between nodes u and v by duv 
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A link is a 3-tuple e = {se, r^,, Pe), where Se € 1^ is the transmitter, re € F is the receiver, and Pe 
is the transmission power In the general setting we allow parallel links with different powers. The set 
of links is denoted by C and m = \C\. We abbreviate and denote the distance dg^r^ by de- Similarly, we 
denote the distance dg^r' by dee'- Note that according to this notation, dee' / c?e'e- 

We use the following radio propagation model. A transmission from point s with power P is received 
at point r with power Pjdl^. The exponent a is called the path loss exponent and is a constant. In most 
practical situations, 2 < a < 6; our algorithm works for any constant a > 0. For links e, e', we use the 
following notation: Se = Pe/d"^ and Se'e — Pe' /d^,^. 

A subset of links L C £ is SiNR-feasible if Se/{N + Y^e'eL\{e} ^e'e) ^ for every e € L. 
This ratio is called the signal-to-noise-interference ratio (SINR), where the constant N is positive and 
models the noise in the system. The threshold /3 is a positive constant. The ratio Se/N is called the 
signal-to-noise ratio (SNR). 

A link e can tolerate an accumulated interference Y^^, Se'e that is at most (^e — jiN) / /3. This amount 
can be considered to be the "interference budget" of e. Let 7e = {l3Se)/{Se — PN). We define three 
measures of how much of the interference budget is "consumed" by a link e'. 

S I 

ae'(e) = ae'(e) = 7e • ae'(e), and de'{e) = min{l, ae'(e)}. 

The value of Og' (e) is called the affectance IIHW09II of the Unk e' on the link e. The affectance is additive, 

so for a set L C C, let aiie) = J2{e'eL:e'jte} ae'(e)- 

Proposition 1 ( ||HW09i ). A set L C £ is smR-feasible iffaiie) < I, for every e e L. 

Following IIHW09II . we define a set L C £ to be ap-signal, if aL{e) < 1/p, for every e £ L. Note 
that L is SiNR-feasible if L is a 1-signal. We also define a set L C £ to be a p-signal, if dL{e) < 1/p, 
for every e € L. Note that L is SiNR-feasible if L is a (1 + e)-signal for some e > 0. 

By Shannon's theorem on the capacity of a link in an additive white Gaussian noise channel IIGal68L 
it follows that the capacity is a function of the SINR. Since we use the same threshold /3 for all the links, 
it follows that the Shannon capacity of a link is either zero (if the SINR is less than /3) or a value 
determined by /? (if the SINR is at least /?). We set the length of a time slot and a packet length so that, 
if interferences are not too large, each link can deliver one packet in one time slot. By setting a unit of 
flow to equal a packet-per-time-slot, all links have unit capacities. We do not assume that /3 > 1; in fact, 
in communications systems /3 may be smaller than one. 

Multi-commodity flows. Recall that a function g : C ^ is a flow from s to t, where s,t G V , 
if it satisfies capacity constraints (i.e., g{e) < 1, for every e G C) and flow conservation constraints in 
every vertex veV\{s,t} (i.e., EeGin(^,) 5(e) = Eeeout(,;) 5(e)). 

We use multi-commodity flows to model multi-hop traffic in a network. The network consists of the 
nodes V and the arcs C, where each arc has a unit capacity. There are k commodities Ri = {si,ii,bi), 
where Si and ii are the source and sink, and bi is the demand of the ith commodity. Consider a vector 
/ = (/i, . . . , fk), where each /, is a flow from Si to tj. We use the following notation: (i) /i(e) 
denotes the flow of the ith flow along e, (ii) equals the amount of flow shipped from Sj to ti, 
(iii) /(e) = Ei=i /i(e)> (iv) |/| = Yli=i \ fil A vector / = (/i, . . . , /fc) is a multi-commodity flow if 
/(e) < 1, for every e G 

We denote by T the polytope of all multi-commodity flows / = (/i, . . . , /fc) such that \ fi\ < hi, 
for every i. For a p > 0, we denote by Tp T the polytope of all multi-commodity flows such that 

m/bi > p. 
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Schedules and multi-commodity flows. We use periodic schedules to support a multi-commodity 
flow using packet routing as follows. We refer to a sequence {Lt}J^Q, where Lt C for each i, as a 
schedule. A schedule is used periodically to determine which links are active in each time slot. Namely, 
time is partitioned into disjoint equal time slots. In time slot t', the links in Lt, for t = t' (mod T) 
are active, namely, they transmit. Each active link transmits one packet of fixed length in a time slot 
(recall that all links have the same unit capacity). The number of time slots T is called the period of the 
schedule. We sometimes represent a schedule S = {Lt}JSQ by a multi-coloring vr : £ ^ 2^'^' - '^^^h 
The set Lt simply equals the preimage of t, namely, Lt = TT~^(t), where vr^^(t) = {e : t £ vr(e)}. 

An SlNR-schedule is a sequence {Lt}JSQ such that Lt is SiNR-feasible for every t. Consider a multi- 
commodity flow / = (/i, . . . , fk) and a schedule S = {Lt}JSQ. We say that the schedule S supports f 
if 

yeeC: T-f{e) < |{t G {0, . . . , T - 1} : e G Lt}| . 

The motivation for this definition is as follows. Consider a store-and-forward packet routing network 
that schedules links according to the schedule S. This network can deliver packets along each link e at 
an average rate of /(e) packets-per-time-slot. 

Buckets and signal diversity. We partition the links into buckets by their received power Se ■ Let 

5min — iiiiiieG£ Se- The ith bucket Bi is defined by 



Bi — |e G /3 I 2 • Smin ^ "^e < 2 ^ • S'min} • 

For a link e G £, define i{e) = [log2(5'e/S'min)J- Then, e G B^^^y The signal diversity o" of £ is the 
number of nonempty buckets. 



Lemma 1. 



a <\a- log2 A + log2 T] 



Proof. Recall that Se = Pe/df^- The signal diversity of C is at most log2(5'max/5'min)> where 
max{5e : e G and S'min = minjS'e : e G C}. Hence, 



log2(5'max/5'min) < log2(%^/^^) 



mm "'max , 

= iog2(r.A"), 

where Pmin = min{Pe : e G /:},Pmax = max{Pe : e G £},dmax = max{de : e G £},dn 
min{(ie : e G C}, as required. □ 



^mm 



Power assignments. In the uniform power assignment, all links transmit with the same power, namely, 
Pf, = Pe' for every two links e and e'. In the linear power assignment, all links receive with the same 
power, namely, Se = Se' for every two links e and e'. 

Assumption on SNR. Our analysis requires that, for every link e, Se/N > (1 + e) • /3, for a constant 
e > 0. Note that if 5e/A^ = /?, then the link cannot tolerate any interference at all, and 7e = oo. 
Our assumption implies that 7e < (1 + e) • /3/e. This assumption can be obtained by increasing the 
transmission power of links whose SNR almost equals f3. Namely, if Se/N /?, then Pe ^ {1 + e) ■ 
Pe. A similar assumption is used in ICKM+OSl . where it is stated in terms of a bi-criteria algorithm. 



Namely, the algorithm uses transmission powers that are greater by a factor of (1 + e) compared to the 
transmission power of the optimal solution. 
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Assumption 1. For every link e £ C, Se/N > {1 + e) ■ /3. 
Proposition 2. Under Assumption\l\ /? < 7e < (1 + e) • /3/e. 

Proof. Recall that7e = g^fgj^ = jiignvTsl- Assumption [T] implies that Se/N > (3. Hence, 7e > (3. 



Se - 1+e- 

/3 



Assumption [T] implies that (3^ < Hence, 



7e 



< 



1 - (3{N/Se) 

(1 + e) ■ (3/e , 



as required. □ 



3 Problem Definition 

The problem MAX THROUGHPUT is formulated as follows. The input consists of: (i) A set of nodes 
V in M? (ii) A set of links C The capacity of each link equals one packet per time-slot, (iii) A set of 
requests {Ri}^^^. Each request is a 3-tuple Ri = {si,ti,bi), where Si £ V is the source, £j € ^ is the 
destination, and bi is the requested packet rate. We assume that every request can be routed, namely, 
there is a path from Sj to ii, for every i € [1..A;]. Since the links have unit capacities, we assume that the 
requested packet rate satisfies hi < n. The output is a multi-commodity flow / = (/i, . . . , fk) G T and 
an SiNR-schedule S = {Lt}JSQ that supports /. The goal is to maximize the total flow |/|. 

The Max-Min Throughput problem has the same input and output. The goal, however, is to 
maximize p, such that / G J^,. Namely, maximize minj=i ^ \ fi\/bi. 



4 Necessary Conditions: siNR-feasibility for links in the same bucket 

In this section we formalize necessary conditions so that a set of links in the same bucket is SINR- 
feasible. In Section [5] we use these conditions to build a LP-relaxation for the problem. 

We begin by expressing ae^(e2) in terms of the distances , dgj, deie2- Note that aei(e2), with 
respect to links that are in the same bucket, depends solely on d^^ and deie2- On the other hand, dei (62), 
with respect to the uniform power model, depends solely on de2 and de^e2- The proof of the following 
proposition is in Appendix |A] 

Proposition 3. 

d, 



62 

62 



V ei , 62 G £ : (62) = I -z 1 in the uniform power model. 

\de1e2 ) 

Throughout this section we assume the following. Let L (1 C denote an SiNR-feasible set of links 
such that all the links in L belong to same bucket Bi. Let e £ Bi denote an arbitrary link (not necessarily 
in L). 

Notation. Define: 

= {e £ L : d^,' < (ie'eji and 
L9 ^{e' eL: d^, > de'e}- 
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4.1 A Geometric Lemma 

The following lemma claims that for every e ^ Bi (not necessarily in L), there exits a set of at most six 
"guards" that "protect" e from interferences by transmitters in L^. 

Lemma 2. There exists a set G of at most six receivers of links in such that 

Ve' e 3c/ G G : d^'g < 2 • d^'e- 

Proof. The set G is found as follows (see Figure [T}: (i) Partition the plane into six sectors centered at 
Te, each with an angle of 60°. Denote these sectors by sector{j), where j G {1, . . . , 6}. (ii) For every 
sector{j), let Cj S denote a link such that the transmitter is closest to among the transmitters 
in sector{j). (iii) Let gj denote a link in such that Vg^ is closest to (If sector{j) lacks transmitters, 
then Qj is not defined). Let G = {vg. denote the set of guards. 



sector(Q) 



sector(5) 



sector (1) 




sector {4) 




sector{2) 



sector{3) 



Figure 1: A depiction of the proof of Lemma |2] 



We first consider the case that e' E is also a guard (e' = gj). In this case choose g = e', and 
d-e'g = de'- But d^' < d^'e sincc e' € L^, as required. We now consider the case that e' £ \ G. Given 
e' € \ G, let j denote the sector that contains Sg'. We claim that d^'g^ < 2 • dg/g. Consider first e' = ej 
(i.e., is a closest sender to rg in sector{j)). Since Vg^ is a closest receiver to Sg^ , we have de^g^ < dg^ . 
Since Cj G L^, we have dg^ < de^e- Thus, (ig^g^ < de^e, as required. 
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Consider now a link e' 7^ ej. The following inequalities iiold: 

de'e > d(..(., (Sej is a closest sender to re) (1) 

de'g, < ds^ise. + dejgj, (triangle ineq. in Ase'Se^Vg^) (2) 

dejgj < de^e, (already proved for ej) (3) 

ds,s,. < de'e- (proved below). (4) 

We now prove Eq. |4](see Figure |2]i. Let s* denote the point along the segment from to Se' such 
that ds'r^ = de e- The triangle Ar^Se-s* is an isosceles triangle. Since Zse r^s* < 60°, it follows 
that the base angle Zr^Se^s* > 60°. Hence, ZreSejSe' > /-r^Se^s* > 60°. Since Zse^VeSe' < 60°, it 
follows that dg^i^se^ ^ de'e, as required. 



sector{l) 




Figure 2: The triangle Ar^Se -s* is an isosceles triangle. The angle r < 60°. The angle 5 > 60°. The 
angle w > 5 > 60°. 

To complete the proof that de'g^ < 2 • dg/g, observe that 

eq.|2] eqs.[3l4] eq.[T] 

□ 

4.2 Necessary Conditions 

Recall that Let L C £ is an SiNR-feasible set of links that belong to same bucket B,;. Let e ^ Bi denote 
an arbitrary link (not necessarily in L). 

Lemma 3. 

^ ae,(e) = 0(l). 
Proof. By Lemma[2l we find a set of "guards" G C L^, such that: 
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(i) IGI < 6, 

(ii) Ve' G 3^ G G : de'g < 2 • d^'e- 
First, let us bound J2e'£L'\G ^eK^)' 



E Me)< E 

dp/ 



<2"+2.^a^,(^), (5) 

where the first line follows from Proposition |3] The second line follows from Lemma |2] The third line, 
again, follows from Proposition |3l 
Since de'{e) < 1, we obtain 

de'ie) < Yl ae'(e) + |G| , (6) 



Hence, 



= Y 7e- de'{e) + \G\ 

< 7e-2"+2.^a^.(5) + |G| 

g&G 

< \G\.[^ + 1 

VmmggGTs 

< 6 ^ ' + 1 



where the first line follows from Equation |6] and the fact that ae'(e) < ae'(e). The second line follows 
from the fact that 7e • ae'(e) = ae'(e). The third line follows from Equation |5l The fourth line follows 
since is SiNR-feasible, that is, ai^i{g) < 1 and dii{g) < l/7g, for every g ^ G. The last line follows 
from Proposition |2] Lemma|2j and |G| < 6. Since, a and e are constants, the lemma follows. □ 

Lemma 4. 

Y aAe) = 0{l). 

e'eLa 

Proof. Pick e* to be a shortest link in L^. It follows from Proposition |3] and the triangle inequality (see 
Figure |3]l that 

Ve' G \ {e*} ■ d /(e*) > - • ( V > - ■ ( — V 

2 \de'e* J ~2 Vde/g + de.g + 4. y 
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Figure 3: The distance dg'e* is depicted by a bold segment. We bound dg'e* by applying the triangle 
inequahty, that is the dashed segments and the length of link e*, de* ■ 

Since e', e* G L^, it follows that de' > de'e and de* > de*e- Since de' > d^ it follows that 

Since 0^9(6*) = 7e* • aL9(e*), it follows: 

Since is SiNR-feasible, it follows that 0^9(6*) < 1. Hence, 

- • — •7e. • ([L^l - 1) < 1^ 

2 3'^ 

\L3\ < 2-3°/7e. +1. 

Proposition [2] implies that = 0(1). Since a is a constant, it follows that \L^\ = 0(1). Since 
Ylie'eL9^e'{e) < |L5|, the lemma follows. □ 

Lemmas |3] and |4] imply the following theorem. 

Theorem 1. Let L denote an SINR-feasible set of links. If L Q Bi, then 

Me&Bi-. ^ aei{e) < aiie) +ae{e) = 0(l). 

{e'£L:d^,>de} 

The following theorem follows from HKeslll Thm 1]. The proof of the following theorem is in Ap- 
pendix |Al 

Theorem 2. Let L denote an SINR-feasible set of links. If L Q Bi, then 

{e'€L:d^,>de} 

5 LP Relaxation 

In this section we formulate the Unear program for the MAX THROUGHPUT and M AX-MiN THROUGHPUT 
problems with arbitrary power assignments. The linear program formulation that we use for computing 
the multi-commodity flow / is as follows. 
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MaxTh^p : F* = maximize |/| subject to 

feT (7) 

VzVeGS, /(e)+ (a^'ic) + ae{e')) ■ f{e') < 1 (8) 

{e'eBi:d^,>de} 



MaxMinTh^^p : R* =maximize p subject to 

/ e -Fp (9) 

Vi Ve G 5, /(e) + ^ (ae'(e) + ae(e')) • fie') < 1 (10) 

Recall that T denotes the polytope of all multi-commodity flows / = (fi,...,fk) such that < 
bi, for every i. Also recall that Tp ^ T for p > denotes the polytope of all multi-commodity flows 
such that \fi\/bi > p. Constraints |71|9] in MaxTh^p and MaxMinTh^p respectively require that the 
/ is a feasible multi-commodity flow with respect to T and Tp. 

Constraints [8l[T0]in MaxTh^p and MaxMinTh^p respectively require that for every bucket 
Bi and for every link e ^ Bi the amount of flow /(e) plus the amount of the weighted symmetric 
interferences is bounded by one. Note that this symmetric interference constraint is with respect to links 
that are longer than e. 

The objective function of MaxTh^p is to maximize the total flow [/[. The objective function of 
MaxMinTh^p is to maximize p, such that / G Tp. Namely, maximize minj=i...fc |/i|/6j. 

We prove on Section |7] that the linear programs MaxTh^p and MaxMinTh^p are relaxations of 
the Max Throughput and Max-Min Throughput problems. 



6 Algorithm 

6.1 Algorithm description 

For simplicity, we assume in this section that all the links are in the same bucket, that is £ C for 
some i. In Section [8] we show how to handle arbitrary power assignment. In Section |9] we extend the 
algorithm so that it assigns limited powers. 

Algorithm overview. We overview the algorithm for the MAX THROUGHPUT problem. Assume for 
simplicity that, C C Bi for some i. First, the optimal flow /* is obtained by solving the linear program 
MAXTHip. We need to find an SiNR-feasible schedule that supports a fraction of /*. Second, we color 
the links using greedy multi-coloring. This coloring induces a preliminary schedule, in which every 
color class is "almost" SiNR-feasible. This preliminary schedule is almost SiNR-feasible since in every 
color class and every link e, the affectance of links that are longer than e on e is at most 1 . However, 
the affectance of shorter links on e may be still unbounded. Finally, we refine this schedule in order to 
obtain an SiNR-feasible schedule. Note that the returned SiNR-feasible schedule supports a fraction of 
the flow /*. We show in Section |7] that this fraction is at least $1(1/ log n). 

Algorithm description. The algorithm for the MAX THROUGHPUT problem proceeds as follows. 

1. Solve the linear program MaxTh^p. Let /* denote the optimal solution. 

2. Remove flow paths that traverse edges with /*(e) < l/(2nm). Let / denote the remaining flow. 
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3. Set T = 2nm. Apply the greedy multi-coloring algorithm greedy-coloring (see Section [631 ) on 
the input ((£, £^), /, d, w, T), where the pair {C, C?) is a complete graph whose set of vertices is 
L, for every link in e G £, d(e) = dg, and w(e, e') = ae(e') + ae'(^) ^ weight function over 
pair of links in L. Let vr : £ ^ 2^'^' ' "^~^J^ denote the computed multi-coloring. 

4. Apply procedure disperse to each color class (7r~^(t)), where t G {0, . . . T — 1}. Let 
denote the dispersed subsets. 

5. Return the schedule {-Lt,i}j=o..T-i,i=i..£(t) and the flow / = (/i,. . . where / = f/{2-£{t)). 

Clearly steps [T]and|5]are polynomial. In Section 163] we show that step|3]is polynomial. In Section [64] 
we show that disperse is polynomial. Therefore, the running time of the algorithm is polynomial. 

Remark 1. The following changes are needed in order to obtain an algorithm for the Max-Min Throughput 
problem: (i) In Item\l}solve the linear program MaxMinTH/^p, (ii) in Item^ remove flow paths that 
traverse edges with /*(e) < 1/ {2'n?km), (Hi) in Item\3]set T = 2v?km. 

6.2 Removing Minuscule Flow Paths 

The greedy multi-coloring algorithm cannot support flows /*(e) < l/(2nm). We mitigate this problem 
simply by peeling off flow paths that traverse edges with a flow smaller than l/(2nm). The formal 
description of this procedure is as follows. (1) Initialize / /. (2) While there exists an edge e with 
/(e) < l/(2nm), remove flow from / until /(e) = 0. This is done by computing flow paths for the 
flow that traverses e, and zeroing the flow along these paths. 

6.3 Greedy Multi- Coloring 

Let G = {V, E) denote an undirected graph with edge weights w : E ^ and node demands 
X : V ^ Assume an ordering of the nodes induced by distinct node lengths d{v). For a set 

V C V, let w{V' , u) = Ylv&v ^(^' "")• Assume that 

Vti G y : x{u) + w{v,u) ■ x{v) < 1 . (11) 

{v&V:d{v)>d{u)} 

Indeed, Constraints [8l [TOl in MaxTh^p and MaxMinTh^p, respectively, imply that the input to 
the greedy coloring algorithm satisfies the assumption in Equation [TT] 

Lemma 5 (Greedy Coloring Lemma). For every integer T, there is multi-coloring vr : F — )• 2^^''''''^~^\ 
such that 

1. VcG{0,...,r-l} VnGvr-Hc):E{.gy:d(.)>d(n)}^K^)<l. 

2. yueV : \tt{u)\ > [x{u)-T\. 

The running time of Algorithm [T] is at most 0{T^ ■ \ V\ ■ \E\). Since T,\E\ and |y| are polynomial, 
it follows that the running time is polynomial. 

Proof. We apply a "first-fit" greedy multi-coloring listed in Algorithm [T] We now prove that this algo- 
rithm succeeds. 
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Algorithm 1 greedy-coloring {{V, E),x, d, w, T) - greedy multi-coloring of V. 



1. Scan the vertices in descending d{v) length order, let u denote the current node. 

(a) C^'i^{cG{0,...,r-l} :u;(^-i(c),u) > 1}. 

(b) If \Cl^'^\ >T -[x{u)- TJ, then return "FAIL". 

(c) 7r(n) ^ first [x{u) ■ T\ colors in {0, . . . , T - 1} \ C^^'^. 

2. Return (tt). 



Let h{u) = \ x{u) ■ T\. Assume, for the sake of contradiction that, \C^'^\ > T — b{u), hence, 

r-6(n) + l < \Cl^'^\ 

< ^ u'(7r-^(c),n) 

< ■ ^{v, u) 

{v:d{u)<d{v)} 

^ b{v)-wiv,u). (12) 

{v:d{u)<d{v)} 

The third line follows from the fact that vertices are scanned in a descending length order, and by a 
rearrangement of the summation order. By adding b{u) to both sides, we obtain: 

T + 1 < [x{u) -Tl + lx{v) -Tl ■w{v,u). (13) 

{v:d{u)<d{v)} 

We divide Eq. [13] by T to obtain a contradiction to Eq. [TT] as required. We conclude, that the greedy 
coloring succeeds, and the lemma follows. □ 

6.4 The dispersion procedure disperse 

The input to the dispersion procedure disperse consists of a set L C £ of links that are assigned the 
same color by the multi-coloring procedure (see Algorithm [T] in Section [63]) . This implies that 

Ve E L : ^ (oe(e') + ae'{e)) < 1. (14) 

{e'eL\{e}:d^,>de} 

The dispersion procedure works in two phases. In the first phase, L is partitioned into 1/3-signal 
sets {Li}i. In the second phase, each subset Lj is further partitioned into 7/6-signal sets {-Lj}^2i- Recall 
that a set of links Lj is SiNR-feasible if Li is a (1 + e)-signal for some e > 0. Since every set in {Lij^^i 
is (7/6)-signal, it follows that every set in {Li}ff[ is SiNR-feasible. 

In Algorithm [2] we list the first phase of the dispersion procedure. Note that if a 1/ 3-signal set J* 
is always found in Line [2a] then L is dispersed into at most log2 \L\ subsets. In Lemma[8]we prove that 
this is indeed possible. 

The second phase follows PHW091 Thm 1]. This phase is implemented by two first-fit bin packing 
procedures. In the first procedure, open 7 bins, scan the links in some order and assign each link to 
the first bin in which its affectance is at most 3/7. In the second procedure, partition each bin into 7 
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Algorithm 2 ^-disperse {L) : partition L Q C into O(logn) 1/3-signaI sets. 

1. i ^ and ^ L. 

2. while L V do 

(a) find a T73-signal set J* C U such that \J'\ > \U\/2. 

(b) L^+i ^ \ J* and i ^ i + 1. 



sub-bins. Scan the links in the reverse order, and again, assign each link to the first bin in which its 
affectance is at most 3/7. 

Proposition |2] implies that step|2]in Algorithmic] terminates after O(logm) iterations. Each of these 
iterations is polynomial. The second phase of the disperse algorithms is clearly polynomial. Therefore, 
the running time of the disperse algorithm is polynomial. 

7 Algorithm Analysis 

In this section we analyze the algorithm presented in Section[6l Recall that it is assumed that all the links 
are in the same bucket, that is £ C Bi for some i. First, we prove that the linear program MaxTh^p 
is a fractional relaxation of the Max Throughput problem. We then show that the greedy coloring 
computes a schedule that supports the flow given by the LP. Unfortunately, this schedule is not an SINR- 
feasible schedule. We then prove that the refinement procedure (Step |4] of the algorithm) generates an 
SlNR-feasible schedule with an O(logn) increase in the approximation ratio. 

Let /* denote an optimal solution of the linear program MaxTh^p, i.e., F* = |/*|. The following 
lemma shows that the linear program MaxTh^p is a relaxation of the Max Throughput problem. 

Lemma 6. There exists a constant A > 1 such that, if S = {Lt}^^ is an SlNR-feasible schedule that 
supports a multi-commodity flow f, then //A is a feasible solution of the linear program MaxTHlp. 
Hence, F* >\f\/\. 

Proof. Clearly //A € J^. Thus, we only need to prove that //A satisfies the constraint in Eq.|8] Consider 
an SlNR-feasible set Lt and an arbitrary link e. By, Theorems [T] and |2j 

{a,,{e) + a,{e'))<0{l). 

{e'eLt:d^,>de} 

It follows that 

E (ae'(e)+ae(e'))<0(l). (15) 

t=0 {e'€Lt:d^,>d^} 

Since /(e') < 7 • |{t : e' G Lt}\, We conclude that 

T-l 

E («e'(e)+ae(e'))> ^ {ae'{e) + a,{e')) ■ f{e'). (16) 

t=0 {e'&Lt:d^,>de} {e'£C:d^,>de} 

Since /(e) < 1, we conclude from Eqs. [Tsl and [T6] that 

/(e) + (ae'(e) + ae(e'))-/(e') <0(1). (17) 

{e'eC:d^,>de} 
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Let A > denote a constant that bounds the left-hand side in Eq. [T7] Then, //A satisfies the constraints 
in Eq. [H as required, and the lemma follows. □ 



Analogously, one could prove also that the linear program MaxMinTh^p is a relaxation of the 
Max-Min Throughput problem. 

Lemma 7. Suppose S = {Lt}J^Q is an SlNR-feasible schedule that supports a multi-commodity flow 
f. If p = minj=i,,jj \ fi\/bi, R* > p/X,for the same constant X > 1 in Lemma^ 

The following proposition gives a lower bound on the optimal throughput. 

Proposition 4. F* > ^ and R* > 

Proof. Without loss of generality, the source and destination of each request are connected. Pick a 
request Ri and a path pi from to tj. Consider the schedule that schedules the links of pi in a round- 
robin fashion. Clearly, this schedule supports a flow / = l/\p\ from Sj to ti along p, where \p\ denotes 
the length of p. This implies that F* > 1/n, as required. The second part of the proposition is proved 
by concatenating k schedules, one schedule per request. The concatenated schedule supports a flow 
f = ifi, ■ ■ ■ , fk), where fi = 1/ (nk) along the path pi. Since bi < n, it follows that \fi\/hi > 1/ {n^k), 
and the proposition follows. □ 

Propositions. |/| > F* /2 

Proof. Let us denote by g the total flow that was removed in step |2] The contribution to the flow 
amount \g\ due to edges with small flow is less than l/(2nm). Since there are m edges, it follows that 
I^I < l/(2?^)■ By Prop. 0] we have F* > i, and the proposition follows. □ 

For the case of MaxMinTh lp, one can show a similar result, that is |/| > R* /2. 

Proposition 6. If T > 2nm, then the greedy multi-coloring algorithm computes a multi-coloring tt 
that induces a schedule that supports f /2. 

Proof. Recall that a schedule S = {Lt}JSQ induced by a multi-coloring vr : £ — 2^'^' ' '^^^J^ is defined 
by Vt : Lt = 7r~^(t), where Tr^^{t) = {e : i G 7r(e)}. Also recall that a schedule S supports / if 
Ve G £ : T ■ /(e) < \{t £ {0, . . . ,T — 1} : e £ Lt}\. Lemma[5]implies that the greedy multi-coloring 
algorithm (see the hsting in Algorithm [U computes multi-coloring vr such that Ve G £ : |vr(e)| > 
[/(e) • TJ. Hence, it suffices to prove that T ■ /(e)/2 < [T ■ /(e)J, for every edge e. Indeed, step|2] 
in the algorithm (see listing in Sec.|6ll implies that if /(e) > 0, then /(e) > 1/T. Let us consider the 
following two cases: (1) If /(e) G [1/T, 2/T), then T ■ /(e)/2 < 1 = [T • /(e)J, (2) if /(e) > 2/T, 
then T ■ /(e)/2 < T ■ (/(e) - 1/T) < [T • /(e)J, as required. □ 

For the case of MaxMinTh lp, one can show the same result if T > 2n?km. 

Lemma 8. If L Q C satisfies Eq. [7?] then there exists a subset J (1 L such that: ( i) J is a 1/ 3-signal, 
and(ii) \ J\ > \L\/2. 

Proof. Define a square matrix A, the rows and columns of which are indexed by L as follows: order 
L in descending length order, so that e' precedes e if d^' > d^. Let A{e,e') = (ae(e') + ae'(e)) and 
A{e, e) = 0. Note that A is symmetric. 

Let A^ denote the upper right triangular submatrix of A. Eq. [14] implies that, 

{e':d^,>de} 
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Hence, the weight of every column in is bounded by 1. This implies that the sum of the entries in 
A'^ is bounded by \L\. By Markov's Inequality, at most half the rows in A'^ have weight greater than 2. 
Let J Q L denote the indexes of the rows in A'^ whose weight is at most 2. Clearly, \ J\ > \L\/2. 

We claim that, for every e G J, the weight of the column is at most 3. Indeed, A(e', e) < 

1. In addition, E{e':d^,<d,} ^(e',e) = E{e':d,, ^<de}^(^'^') — 2 since this is the sum of the row indexed 
e in A^. This implies that aj{e) < 3, for every e E J, and the lemma follows. □ 

Proposition 7, The dispersion procedure partitions every color class vr^^(t) into O(logm) SINR- 
feasible sets. 

Proof. Recall that the dispersion procedure disperse consists of two phases. The first phase is the i- 
disperse{iT~^{t)) algorithm (see the Usting in Algorithm IS, and the second phase is implemented by 
two first-fit packing procedures. 

Let us consider the first phase. NotethatL° = TT^^{t). Since < 1^1/2,1^1^11 ^-disperse{Tr~^{t)) 

requires at most log2 |7r^^(t)| iterations. Hence, it partitions vr^^(t) into at most log2 |7r~^(t)| sets, 
where each set is a 1/ 3-signal set. 

Now, in the second phase each of these sets is partitioned into 49 subsets. The lemma follows. □ 

Theorem 3. IfAssumption\l}holds, and all the links are in the same bucket, then there exists an 0(log n)- 
approximation algorithm for the MAX THROUGHPUT and the MAX-Min THROUGHPUT problems. 

Proof. Let OPT denote the maximum total throughput. By Lemma [6l F* > opt/ A = O(opt). Recall 
that /* denotes an optimal solution of MaxTHlp. By Prop. [5] [/| > |/* |/2, and by Prop. |6l the multi- 
coloring TT supports //2. By Prop. |7J the dispersion procedure reduces the throughput by a factor 
of O(logm). Since there are no parallel edges, logm = O(logn). Thus, the final throughput is 
|/|/0(logn) = OPT/0(logn), and the theorem follows. □ 

Since in the Unear power assignment all links receive with same power, all the links are in the same 
bucket. We conclude with the following result for the linear power assignment. 

Corollary 4. If Assumption \J} holds, then there exists an O (log n)-approximation algorithm for the 
Max Throughput and the Max-Min Throughput problems in the linear power assignment. 



8 Given Arbitrary Transmission Powers 

In this section we show how to apply the algorithm presented in Section|6]to the case in which transmis- 
sion power Pe of each link e is part of the input. Note that Pe may be arbitrary. 

Theorem 5. If Assumption\l\holds, then there exists an 0(log n ■ (log A + log T)) -approximation algo- 
rithm for the Max Throughput and the Max-Min Throughput problems when the link transmis- 
sion powers are part of the input. 

Proof sketch: We construct an SiNR-feasible schedule and its supported flow. The construction pro- 
ceeds as follows: (1) solve the matching LP, (2) remove the minuscule flow paths as described in Item|2j 
(3) run Items [3][5] for every bucket separately, (4) concatenate the output schedules, to obtain an SINR- 
feasible schedule of all the links in £. Step (3) of this construction reduces the flow by a factor of at most 
0(log n). Step (4) of this construction reduces the flow by an additional factor of at most the number of 
nonempty buckets, that is 0(log A + log F). □ 
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9 Limited Powers 



In this section we consider the case in which the algorithm needs to assign a power to each link. 
The assigned powers must satisfy Pmin < ^'e < ^max- To simplify the description, assume that 
log2(Pmax/-Pmin) IS an integer, denoted by i. 

We reduce this problem to the case of given arbitrary powers as follows. For each pair of nodes 
(n, v), define £ + 1 parallel links, where the transmission power of the ith copy equals 2* • Pmin- 

Theorem 6. Assume that, for every link e, (Pmin/ d^) /N > (1 + e) • /3. Then, there exists an 0((log n + 

log log r)-(log A+log T)) -approximation algorithm for the MAX THROUGHPUT and the Max-Min THROUGHPUT 

problems when the link transmission powers are in the range [Pmin, ^max]- 

Proof sketch: Note that the number of links increases by a factor of O(logr). This implies that the 
log n factor increases to (log n + log log F). 

The important observation is that there exists a solution that uses the discrete power assignments 
2* • Pe and achieves a throughput that is a constant fraction of the optimal throughput. The theorem 
follows then from Theorem |5] 

The proof of this observation proceeds as follows. Given an optimal schedule, refine each time slot 
so that it is a p-signal for p = 2. This reduces the throughput only by a constant factor (see IIHW09[ 
Thm 1]). Round up each transmission power to the smallest discrete power that satisfies Assumption [1] 
This increases the affectance by at most a factor of two, thus the resulting schedule is SiNR-feasible. 
Moreover, the schedule uses links with powers that satisfy Assumption [T] □ 
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A Proofs 

Proposition |3j 



^ \ '-''6162 / \ ""61 



162 



V ei, 62 E £ : (62) = ( ) in the uniform power model. 

,"^6162 , 
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Proof. Recall that cieK^) — '^e — Pe/d", and 5e'e = -Pe'/'^e'e- Note that every two hnks ei, 62 G 
Bi, satisfy that Se^/Se^ G (1/2, 2). Hence, 

aei(e2) = 



Seie2 _ 


Seie2 


Se2 








Se2 




a Q 


\de1e2 J 


5*62 



as required. 

On the other hand, in the uniform power model assignment, all links transmit with the same power, 
namely Pg = Pe' for every two links e and e'. Hence, 



aei(e2) 



Seie2 

Se2 

-Pei/<^eie2 

Pe2 /de2 
1 \a 
"62 1 



de\e2 



as required. □ 
Theorem |2] Let L denote an SINR-feasible set of links. IfLCl Bi, then 

yeGBi-. Yl = Oil)- 

{e'€L:d^,>de} 

Proof. Theorem 1 in HKeslU implies that 



It follows that. 



{e'€L:d^,>d,} ^ \ ce / J 

> mmll,y ae{e')\ 

{e'€L:d^,>de} ^ ^ 

{e'eL:d^,>de} ^ 

^ "^^"{''2.(l + .)./3-"^(^')} 



{e'&L:d^,>de} 

where the second line follows since L Q Bi and Proposition [3] The third line follows from the definition 
of ae(e')- The last line follows from Proposition |2] The theorem follows, since 2.(i+e)./3 = 0{1) and 
since ae'(e) = min{l, ae'(e)}. □ 
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